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Abstract 

We study the longitudinal motion of beam particles under the action of a single 
resonator wave induced by the beam itself. Based on the method of multiple 
scales we derive a system of coupled amplitude equations for the slowly varying 
part of the longitudinal distribution function and for the resonator wave envelope, 
corresponding to an arbitrary wave number. The equation governing the slow 
evolution of the voltage envelope is shown to be of Ginzburg-Landau type. 
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1 Introduction 



So far, extensive work has been performed on the linear stability analysis of collective motion 
in particle accelerators [|T|. Nonlinear theories 0—0 of wave interaction and formation of 
patterns and coherent structures in intense beams are however less prevalent, in part, due 
to the mathematical complexity of the subject, but also because of the commonly spread 
opinion that highly nonlinear regime is associated with poor machine performance that is 
best to be avoided. 

Nevertheless, nonlinear wave interaction is a well observed phenomenon 0, |] in present 
machines, complete and self-consistent theory explaining the processes, leading to the forma- 
tion of self-organized structures mentioned above is far from being established. The present 
paper is aimed as an attempt in this direction. 

The problem addressed here (perhaps, the simplest one) is the evolution of a beam in 
longitudinal direction under the influence of a resonator voltage induced by the beam itself. 
Linear theory is obviously unable to explain bunch (droplet) formation and bunch breakoff 
(especially in the highly damped regime), phenomena that have been observed by numerical 
simulations @, @, 0, but it should be considered as the first important step towards our 
final goal - nonlinear model of wave interaction developed in Section 3. 

It is well-known that within the framework of linear stability analysis the solution of the 
original problem is represented as a superposition of plane waves with constant amplitudes, 
while the phases are determined by the spectrum of solutions to the dispersion equation. 
Moreover, the wave amplitudes are completely arbitrary and independent of the spatial and 
temporal variables. The effect of nonlinearities is to cause variation in the amplitudes in 
both space and time. We are interested in describing these variations, since they govern the 
relatively slow process of formation of self-organized patterns and coherent structures. 

The importance of the linear theory is embedded in the dispersion relation and the type 
of solutions it possesses. If the dispersion relation has no imaginary parts (no dissipation of 
energy occurs and no pumping from external energy sources is available) and its solutions, 
that is the wave frequency as a function of the wave number are all real, then the correspond- 
ing amplitude equations describing the evolution of the wave envelopes will be of nonlinear 
Schrodinger type. Another possibility arises for conservative systems when some of the roots 
of the dispersion equation appear in complex conjugate pairs. Then the amplitude equations 
can be shown to be of the so called AB-type |TD[ . For open systems (like the system studied 
here) the dispersion relation is in general a complex valued function of the wave frequency 
and wave number and therefore its solutions will be complex. It can be shown ]10[ that 
the equation governing the slow evolution of the wave amplitudes in this case will be the 
Ginzburg-Landau equation. 

Based on the renormalization group approach we have recently derived a Ginzburg- 
Landau equation for the amplitude of the resonator voltage in the case of a coasting beam 
0. The derivation has been carried out under the assumption that the spatial evolution of 
the system is much slower compared to the temporal one. This restriction has been removed 
here, and the present paper may be considered as an extension of |J. 

Using the method of multiple scales we derive a set of coupled amplitude equations for 
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the slowly varying part of the longitudinal distribution function and for the intensity of a 
single resonator wave with an arbitrary wave number (and wave frequency, specified as a 
solution to the linear dispersion equation). The equation governing the evolution of the 
voltage envelope is shown to be of Ginzburg-Landau type. 



2 Formulation of the Problem 

It is well-known that the longitudinal dynamics of an individual beam particle is governed 
by the set of equations || 



dzi 7 a T7i dAE eu s V RF . eu s 

— = k Q AE ; — — = — (srn0-sm0 s ) + — V 1: (2.1, 

dt dt 2ir 2tc 



where 



is the proportionality constant between the frequency deviation of a non synchronous particle 
with respect to the frequency uo s of the synchronous one, and the energy deviation AE = 
E — E s . The quantity k also involves the phase slip coefficient r] = cum — 7<T 2 , where % is 
the momentum compaction factor M. The variables 



Z\ = 9 — uj s t ; <fi — <f) s — hz\. (2.3) 

are the azimuthal displacement of the particle with respect to the synchronous one, and the 
phase of the RF field, respectively. Here Vrf is the amplitude of the RF voltage and h is 
the harmonic number. Apart from the RF field we assume that beam motion is influenced 
by a resonator voltage V\ due to a broad band impedance 



92Vl 2 9Vl I ^y^^eTZdh 
dz\ dz\ u s dt ' 



where 



uj = — ; 7 = 7^ ; hiO',t)= fdAE{u s + k AE)f 1 {e,AE;t), (2.5) 
lo s 2Q J 

fi(9, AE; t) is the longitudinal distribution function, uj t is the resonant frequency, Q is the 
quality factor of the resonator and 1Z is the resonator shunt impedance. 

It is convenient to pass to a new independent variable ( "time" ) 9 and to the new dimen- 
sionless variables @, @: 
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r = u s 9 ; z = Zi^M ; u = —— , (2.6) 

/v. u. 



f 1 (9,AE;t) = ^Lf(z,u;9) ; V X = \ X V ; h=ou s p I, (2.7) 



where 



v s = ; Ai = 2>yQeR,uj s pQ. (2.8) 

In the above expressions the quantity po is the uniform beam density in the thermodynamic 
limit. The linearized equations of motion ( |2.1|) and equation (|2.4j ) in these variables read as: 

dz du , 

— = u ; — = -z + AF, (2.9) 
ar ar 



<9 2 y <9V 9/ /• 

-2 7o — + U *V=-— ; I(z;6) = J du(l + Uy/^)f(z,u;9), (2.10) 



dz 2 
where 



7 w , e 2 n-f k p 
7o = -p= ; = -= ; A = — -. (2.11 

IV. HVc^Ve 



We can now write the Vlasov equation for the longitudinal distribution function f(z,u;9), 
which combined with the equation for the resonator voltage V(z; 9) 

^- + «« ^TT" + = °' ( 2 - 12 ) 

ar 02; au au 



I(z;9) = I du(l + u^)f(z,u;9), (2.14) 
comprises the starting point for our subsequent analysis. 
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3 Derivation of the Amplitude Equations for a Coast- 
ing Beam 

In this Section we analyze the simplest case of a coasting beam. The model equations ( |2.12|) 
and fl2.13| ) acquire the form H 



2 7ir- + ^ V= / d« — -— , (3.2) 



<9z 2 J \89 dz f 

where the parameter A should be calculated for v s — 1. In what follows it will be convenient 
to write the above equations more compactly as: 



where we have introduced the linear operators 



c(°- (37) 



(Q(z,u;6)} = I dug(z,u;6). (3.8) 
To obtain the desired amplitude equation for nonlinear waves we use the method of multiple 



scales [[10[| , [II] . The key point of this approach is to introduce slow temporal as well as 



spatial scales according to the relations: 
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6 ; T x = e9 ; T 2 = e 2 ; ... ; T n = e n 6 ; ... (3.9) 



z ; zi = ez ; z 2 = e z ; ... ; z n = e n z ; ... (3.10) 

where e is a formal small parameter. Next is to utilize the perturbation expansion of the 
longitudinal distribution function /, the resonator voltage V 



DC 



f = fo(u) + J2e k f k ■ V = J2e k V k , (3.11) 

k=l k=l 



and the operator expansions 



<9T fc 9z ^ «^ fc y 



where V 2 implies differentiation with respect to d/dz. Substituting them back into (|3.3| ) 
and (|3.4|) we obtain the corresponding perturbation equations order by order. It is worth 
noting that without loss of generality we can miss out the spatial scale z 2 , because it can be 
transformed away by a simple change of the reference frame. For the sake of saving space 
we will omit the explicit substitution and subsequent calculations and state the final result 
order by order. 

First order 0(e): 

Ffi+XV 1 ^ = 0, (3.15) 
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W = C(h). 

Second order 0(e 2 ): 



(3.16) 



Th + W 2 ^ = -hh ~ AVi^, (3.17) 



VV 2 = C(f 2 ) + A(/x) - V^. (3.18) 



Third order 0(e 3 ): 



?h + AV 3 ^ = -T x h - hh -wM^- \vS±, (3.19) 

OU OU OU 



VV 3 = C(h) + Ci{f 2 ) + A(/i> - V z -^- - (3.20) 



where T n and £ n are the corresponding operators, calculated for T n and z n . 

In order to solve consistently the perturbation equations for each order we need a unique 
equation for one of the unknowns; it is more convenient to have a sole equation for the 
distribution functions /„ alone. This will prove later to be very efficient for the removal 
of secular terms that appear in higher orders. By inspecting the above equations order by 
order one can catch their general form: 

4 + AK^ = «n ; VV n = C(f n ) + p n , (3.21) 

where a n and (3 n are known functions, determined from previous orders. Eliminating V n we 
obtain: 

^ fn + A |!r^ /ri) = ~ X ^ Pn + ^ (3 ' 22) 

Let us now proceed with solving the perturbation equations. The analysis of the first order 
equations (linearized equations) is quite standard, and for the one-wave solution we readily 
obtain: 



V 1 = E{z n -T n y* + E*{z n -T n )e-^\ (3.23) 
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/i = -A 



dfo 
du 



E(z n ; T n ) 
!F(iQ, —ik, u) 



+ 



E*(z n ;T n ) 



-up 



-ik, u) 



+ F(z n ,u;T n ) 



(3.24) 



with 



<p = QO- kz, (3.25) 

where given the wave number k, the wave frequency fi(fc) is a solution to the dispersion 
equation: 



V(k, Q(k)) = 0. (3.26) 

The dispersion function D{k, Q) is proportional to the dielectric permittivity of the beam 
and is given by the expression 

V{k, O) = V(-ik) + A£(in, -<*)( _J— (3.27) 

where 



Fe i<p = F{itt,-ik,u)e i<p ; Ve** = V(-iJfe)e^ ; = -ifc)e iv . (3.28) 

Note that the wave frequency has the following symmetry property: 

Q*(k) = -Q(-k). (3.29) 



The functions E(z n ;T n ) and F(z n ,u;T n ) in equations (|3.23|) and (|3.24|) are the amplitude 



function we wish to determine. Clearly, these functions are constants with respect to the 
fast scales, but to this end they are allowed to be generic functions of the slow ones. 

In order to specify the dependence of the amplitude functions on the slow scales, that 
is to derive the desired amplitude equations one need to go beyond the first order. The 
first step is to evaluate the right hand side of equation ( |3.22| ) corresponding to the second 
order with the already found solution ( |3.23j ) and ( [3.24 ) for the first order. This yields terms 



(proportional to e lip ) belonging to the kernel of the linear operator on the left hand side 
of equation (|3.22| ), which consequently give rise to the so called secular contributions to 
the perturbative solution. If the spectrum of solutions to the dispersion equation ( |3.26| ) is 
complex (as is in our case), terms proportional to e - 2Im ( n ) e appear on the right hand side 
of ( |3.22|) . Since, the imaginary part of the wave frequency we consider small, the factor 
e -2/m(«)0 j g s i ow iy varying in 9 and we can replace it by e - 2Im ( n ) T ™ ( where the slow temporal 
scale T n is to be specified later. This in turn produces additional secular terms, which need 
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to be taken care of as well. (Note that exactly for this purpose we have chosen two amplitude 
functions at first order). The procedure to avoid secular terms is to impose certain conditions 
on the amplitudes E(z n ;T n ) and F(z n ,u;T n ), that guarantee exact cancellation of all terms 
proportional to e tip and terms constant in the fast scales z and 6 (containing e 



-2Im(n)T n 



on 



the right hand side of equation ( |3.22j ). One can easily check by direct calculation that the 
above mentioned conditions read as: 



dV dE dV dE 



dVt dT x dk dzi 



-iXC 




(3.30) 



f x F + 2X 2 Im(Q) 



9 1 df Q \. pl2 



du 



T 



du 



£\2 e -2Im(n)T n 



±dfo 
uj 2 du 



C X {F). 



(3.31) 



Noting that the group velocity of the wave Q g = dQ/dk is given by 



dV dV dtt _ 
~dk + dfllk~ 



a. 



dV (dV^ 

~^k[dn 



(3.32) 



we get 



dE 



-Ki) 



(3.33) 



The above equations ( |3 . 3 1| ) and ( |3.33| ) are the amplitude equations to first order. Note 
that if Im(Q) = we could simply set F equal to zero and then equation ( |3.33| ) would 
describe the symmetry properties of the original system (|3.1|) and ( p.2|) with respect to a 
linear plane wave solution. However, we are interested in the nonlinear interaction between 
waves (of increasing harmonicity) generated order by order, and as it can be easily seen the 
first nontrivial result taking into account this interaction will come out at third order. To 
pursue this we need the explicit (non secular) second order solutions for / 2 and V-i- 

Solving the second order equation (|3.22j) with the remaining non secular part of the 



second order right hand side and then solving equation ( 3.18 ) with the already determined 
fi we find 



h = S F (k, Q, u)E 2 e 2iip + c.c. + F 2 (z n , u; T n 



(3.34) 



V 2 = S v (k, n)E 2 e 2i * + f v e^ + c.c. + G v (z n , T„; [F]), 



(3.35) 
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where c.c. denotes complex conjugation. Without loss of generality we can set the generic 
function F 2 (z n , u; T n ) equal to zero. Note that, in case Im(Q) = we could have set F = 0, 
as mentioned earlier, but we should keep the function F 2 nonzero in order to cancel third 
order secular terms depending on the slow scales only. Moreover, the functions Sf, Sv, fv 
and the functional Gy([F]) of the amplitude F are given by the following expressions: 



X 2 V{-2ik) 







2 V{2k,2VL)F{iVL,-ik,u) du 



dfo 



J 7 (iQ, -ik, u) du 



(3.36) 



S v (k,Q) = X 



2 C{iVL, —ik) I 1 







V(2k,2Vt) \F(iQ,-ik,u)du 



dfo 



T(iVL, -ik,u) du _ 



(3.37) 



V(-ik) 



iXXLxE - V k (-ik) 



dE 

dzi 



(3.38) 



G v (z n ,T n ;[F]) = —£ 1 (F), 



(3.39) 



i(k,n) 



dfo' 



^(iQ, -ik,u) du /' 



(3.40) 



where the fc-index implies differentiation with respect to k. 

The last step consists in evaluating the right hand side of equation (|3.22 ), corresponding 
to the third order with the already found first and second order solutions. Removal of 
secular terms in the slow scales leads us finally to the amplitude equation for the function 

F(z n , it; T n ), that is 



d f 2 dfo\ 2A7<9/ d p dF ~ 



X 2 u 2 



d (± d AfvE* + l(l d 4V V E 



du \F* du 



du yj 7 du 



-2Im(tt)T 2 



(3.41) 



Elimination of secular terms in the fast scales leads us to a generalized cubic Ginzburg- 
Landau type of equation for the amplitude E(z n ,T n ): 



X 2 CG v {[F])E + XCG{[F})f 



v, 



(3.42) 
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where the coefficients a(k), A(k), B(k) and C(k) are given by the expressions: 



«(§) 



(3.43) 



A(k) 



1 + -J- 
V 



V k + i\X{l + tt 6 



(3.44) 



B(k) = C{ 



1 dS F 
T du 




(3.45) 



C(k) 



\Fdu\F du 



and the functional {7([-F]) of the amplitude F can be written as 



(3.46) 



<3{[F]) 




(3.47) 



Equations ( [3.41|) and ( |3.42|) comprise the system of coupled amplitude equations for 
the intensity of a resonator wave with a wave number k and the slowly varying part of 
the longitudinal distribution function. Note that the dependence on the temporal scale 
Ti (involving derivatives with respect to 7\) in equations ( pj.41| ) and ( |3.42| ) through the 
operator L\ and the function fy can be eliminated in principle by using the first order 
equations ( |3.31| ) and ( |3.33j ). As a result one obtains a system of coupled second order partial 
differential equations for F and E with respect to the variables T2 and Z\. 



4 Concluding Remarks 

We have studied the longitudinal dynamics of particles moving in an accelerator under the 
action of a collective force due to a resonator voltage. For a sufficiently high beam density 
(relatively large value of the parameter A) the nonlinear wave coupling, described by the 
nonlinear term in the Vlasov equation becomes important, and has to be taken into account. 
This is manifested in a spatio-temporal modulation of the wave amplitudes in unison with the 
slow process of particle redistribution. As a result of this wave-particle interaction (coupling 
between resonator waves and particle distribution modes) coherent, self-organized patterns 
can be formed in a wide range of relevant parameters. 

We have analyzed the slow evolution of the amplitude of a single resonator wave with an 
arbitrary wave number k (and wave frequency fl(k) defined as a solution to the dispersion 
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relation). Using the method of multiple scales a system of coupled amplitude equations for 
the resonator wave envelope and for the slowly varying part of the longitudinal distribu- 
tion function has been derived. As expected, the equation for the resonator wave envelope 
is a generalized cubic Ginzburg-Landau (GCGE) equation. We argue that these ampli- 
tude equations govern the (relatively) slow process of formation of coherent structures and 
establishment of wave-particle equilibrium. 
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